Fate of Yang-Mills black hole in early Universe 
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Abstract. According to the Big Bang Theory as we go back in time the Universe becomes progressively hotter and denser. 
This leads us to believe that the early Universe was filled with hot plasma of elementary particles. Among many questions 
concerning this phase of history of the Universe there are questions of existence and fate of magnetic monopoles and pri- 
mordial black holes. Static solution of Einstein- Yang-Mills system may be used as a toy model for such a black hole. Using 
methods of field theory we will show that its existence and regularity depend crucially on the presence of fermions around it. 
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INTRODUCTION 

Recent observational data strongly suggest the existence 
of dark matter (DM) 1 1 ] . One of a few remaining DM 
candidates compatible with the Standard Model (SM) 
are primordial black holes (PBH) |2]. These black holes 
were formed in the early Universe during gravitational 
collapse of density fluctuation. The existence of the PBH 
was first proposed in |i3i, |4|, iSJ and their formation dur- 
ing inflation or phase transition in the early Universe was 
discussed in iSSSl- Another type of objects that could 
be formed during phase transition in early Universe are 
topological defects like cosmic strings and magnetic 
monopoles ["9^. If these objects really formed during evo- 
lution of the Universe we should consider the possibil- 
ity of interactions among them. Within the framework 
of field theory a simple model of the system that allows 
the existence of magnetic monopoles may be given by 
Yang-Mills (YM) theory with SU{2) gauge group. Parti- 
cle like solutions of Einstein- Yang-Mills (EYM) system 
describing the magnetic monopole were first constructed 
in ifToi] . On the other hand, solutions describing a black 
hole with magnetic monopole in this context were pre- 
sented in ifTTIl . There was also shown in lfT2ll that finite- 
ness of the black hole mass demands absence of an elec- 
tric part of YM field. 

Another important component of the early Universe 
whose interaction with PBH we should consider is 
fermionic matter. To describe these interactions we use 
methods of field theory. As a toy model of PBH we take 
the aforementioned spherically symmetric YM black 
hole. As a representation of fermionic matter we use 
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Dirac field. Using dimensional reduction and bosoniza- 
tion we will analyze the influence of fermions on the con- 
sidered PBH. 



DIRAC EQUATION IN MAGNETIC 
YANG-MILLS BLACK HOLE 
BACKGROUND 

A general form of a spherically symmetric and time 
independent metric is 

ds^ = ~A^{r)dt^ + B-^{r)dr^ + r^dQ.^ , (1) 

where r is an usual radial coordinate and dQ.^ — dO^ + 
sin^(0)(i(/)^ is a standard metric on two dimensional 
sphere. The location of an event horizon of a black hole 
is given by the largest positive root of A^{r), which is 
denoted as r//. Actually for analyzing fermions equations 
of motion it is more convenient to introduce the so-called 
tortoise coordinate r*. This coordinate is given by the 
following relation: 

B(r) 

dr, = -^dr (2) 
A{r) 

Our metric expressed in coordinates (f , , , 0) is 

ds^ =A^{-dt^ + drl)+C^dQ.^, (3) 

where we introduce a function C to keep in mind that 
now we have r ~ '"(?"*)• From now on, to simplify our 
notation we will omit explicit writing of arguments of 
functions. During our computations we will use tetrad 
formalism. The local ortogonal tetrad is defined by the 
relation = ^'j^^jitlij^ where e'^ is an element of the 
tetrad, g^y and 7],^ are metric tensors in curved and flat 



spacetimes respectively. We use the following signature 
convention: r;oo = -1, J?ii = T722 = 1733 = +1- 

The general form of static and spherically symmetric 
Yang-Mills potential may be written as 



1 - w(r) 
IXC 



(4) 



where a, = (aq , a 1 ) and w represent electric and magnetic 
parts of YM field, «' is a unit vector normal to the sphere, 
A is a coupling constant, Ejjk stands for a totally anti- 
symmetric Levi-Civita symbol, and T*^ are generators of 
SU (2) gauge group represented by Pauli matrices. As the 
representation for four dimensional gamma matrices we 
choose Weyl basis 
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-C7' 



a' 




(5) 



where / is a unit two dimensional matrix and cr' are three 
Pauli matrices. Gamma matrices in curved spacetime are 
given by the standard relation: y'' = e'^y'^. Representing 

tion as 



a spinor Xj/ as Xj/ = 



we may write the Dirac equa- 



ip' \(fL - m\(/R = 0, (6) 
where operators p are given 



w— 1 



p =A dt-iX{[aao±(yai]n-X±——n-GXT}+ 

zaC 



± a ■ nA- ' dr, ± a ■ n{A-'C-'dr,C+ i^^-'A' ' dr,A)~ 

(7) 



±C-^Dsi. 



In the above formula a bar over a quantity represents 
three dimensional vector, a dot - scalar multiplication, x 
- vector multiplication, and D^i is Dirac operator on S^. 
We are mainly interested in fermions in s-wave sector, 
which is described by the lowest eigenvalue of D^i . The 
action of operators h ■ fj_n_- ct x t, a • « and D^i on these 
states may be found in II13I1 . We may use this knowledge 
to integrate over the angles in four dimensional Dirac 
action to obtain effective two dimensional theory in (f , r*) 
plane. 

On the other hand, dimensional reduction of the YM 
action gives us the following Lagrangian: 



C2 



LyM-2c, = -^.fabf" Mw^l - ^(|W|^ - 1)^ 



1 



(8) 



where da = V„ - f^t = daBb - dtB,,, and B„ = e[,ai. 
From now on latin letters from the beginning of alphabet 
will label indexes from curved two dimensional (f,r*) 
spacetime and those from the middle from flat spacetime. 



MASSLESS FERMIONS 

In the massles case we see from (|6]l that chiralities decou- 
ple and effective two dimensional theory contains two 
fermionic fields connected to y/i and xifg. The effective 
Lagrangian for the field is 

^FR-2d - -iCRfVaGR - XB,GRf'fGR+ 

+VGRfLGR-VGRfRGR, (9) 

where V = ^. In this formula Gr is two dimensional 
spinor field connected to Yr by rescaling and multipli- 
cation by appropriate a matrices ifisll . Two dimensional 
flat spacetime gamma matrices are given by the follow- 
ing relations: 



{f ,f } = 2tj'^', rjoo = -l = -T?ii, 

f = ff = (y\ ?L/R = Ui±f)- (10) 



The Lagrangian for the field connected to differs from 
(|9|l only by a sign of a term proportional to Ba ■ 

Analyzing equations of motion for fermions in curved 
spacetime is a highly nontrivial and difficult task. How- 
ever, in case at hand since we can express our prob- 
lem in form of effective two dimensional theory we 
may use bosonization technique |16]. This technique is 
well defined for flat spacetime two dimensional prob- 
lems and also for asymptotically flat spacetimes ifTTlfTsll . 
The merit of bosonization is that we express a fermionic 
sector of our theory in terms of complex scalar field. Ba- 
sic bosonization formulas in our case are given below 



(11) 



(12) 



\l/YR\l/^be-^''^^. (13) 



After bosonization from (|9) we obtain the following 
scalar Lagrangian: 



^BR = -^^a^R^^R - XBa^y^R- 

2 J% 



^Vb{e 



), (14) 



from which we derive an equation of motion for ^r field 



V„V"0K + ^V„B"+4//7V7rycos(2V^0fi) =0. (15) 

Then, we derive an equation of motion for scalar field 
connected to the i/z, sector of our original fermionic 



theory analogically: 
/L 

V„V«0i- ^V«B" + 4^•/7V7rycos(2^/7r0z.) =0. (16) 

Equations ( fTSl l and ( fT6b are highly nonlinear and we 
were not able to find a solution in terms of known special 
functions in the whole spacetime. Nevertheless, we may 
find some useful information about their solutions by 
analyzing them in two asymptotic regions, namely in 
near horizon region and in large region. 

In the large region we have that 

A^wl, C = r, w«±l, r*-r. (17) 
Taking this into account the equation (fTSl l turns into 

-d^(j)R + dl(j)R + 4iby/7i^-^COs{2y/7Z(j)R)^0. (18) 

After dropping a term proportional to (9(r^') we get 
a free wave equation 

-d^<l)R + dl(j)R=0. (19) 

A regular solution to this equation may be expressed as 
a plane wave 

fe=coe-""('±^'), (20) 

where cq is an integration constant. On the other hand, 
in the near horizon region we have 

A^^2K{r-rH), C^kh. r - r,, ^ e^"'' , (21) 

where k is surface gravity of a Yang-Mills black hole. 
The equation ( fTsT i in this region takes the form 

-d^<l>R + ^R + i^t^^^^,2.r, eos(2y^,?)«) = 0. 

(22) 

Because as we approach a black hole horizon — oo, 
we may drop a term proportional to e^^''* in the above 
equation. From this we see that a regular solution is again 
a time dependent plane wave given by 

(/)fi-cie-"»('±^'), (23) 

where ci is some other integration constant. On the basis 
of this analysis we may conclude that a regular solution 
to the considered scalar field equation will be time de- 
pendent. But both scalar fields represent fermionic cur- 
rents and their time dependence ultimately means that 
fermionic fields will also be time dependent. To see how 
this may influence YM field we use equations of motion 



for electric and magnetic parts of this field in the pres- 
ence of fermions. After using bosonization formulas (ITTl 
these equations read iTsIl 

V„ [C'-f] - 2 \w\^b' = A [vV« - V^0i] , (24) 

V„V«w~ ^w(|w|2 - l)+2wBaB" = 
2b 

i— [sin(2y7r0ff) + sin(2\/7r0L)]. (25) 

Conclusions that stem from the above equations are 
the following. First, let us remind that equations (fTST i 
and ( fT6b differ only by sign in front of the term propor- 
tional to Ba and we assume that initially a black hole 
has only magnetic charge (B„ = 0). Second, from equa- 
tion (l24l i we see that in this case contributions of scalar 
fields cancel each other and B„ = is still a valid so- 
lution to ( l24l i. Third, from equation (|25T l we see that 
there is a nonzero contribution of scalar fields to mag- 
netic part of YM field. But since these fields are time de- 
pendent so should be the resulting w. On the other hand, 
through Einstein equations, this results in time depen- 
dence of metric tensor elements and ultimately means 
that the assumption of staticity of Yang-Mills black hole 
is destroyed in the presence of massless fermions. 

MASSIVE FERMIONS 

For massive fermions we see from (|6) that chiralities are 
mixed up. To use bosonization in this case we need to 
make an additional assumption about the form of y/^ 
and We use the following ansatz: 

Gl = iO^Gr = G, (26) 

where Gl and are two dimensional spinor fields con- 
nected with y/L and yfR respectively |15]. Having this in 
mind an effective two dimensional fermionic Lagrangian 
is as follows: 

■^GF-2d = -iGfVaG - XBaGffG^ 

+ (y + m)G7LG+(m-y)G7RG, (27) 

where like in massless case V = ^. Using the same 
bosonization formulas as in the massless case we arrive 
at the following scalar Lagrangian: 

2 ^Jit 
+ (y + m)he-'-^^ + (m - y)be-'^'-^^ . (28) 



The equation of motion for scalar field in this case is 

+ 4ji7\/7r|ycos(2y7r0) + /msin(2y7r0)| =0. (29) 

In the large limit, after dropping terms proportional to 
(9(r^'), we obtain a sine-Gordon equation 

-d}(^ + dl <j) - Ab^/nm sin(2y7r0 ) = 0, (30) 

for which a regular time dependent and decaying at 
infinity (kink type) solution is given by 

2 / - / SfcT'/' fr -ufA 

= -^arctan|^e VT^^ * M. (31) 

On the other hand, the same type of analysis like in 
the massless case revealed that in the near horizon limit 
we again have a free wave equation with a regular time 
dependent solution in form of a plane wave: 

=C2e"''"('±'-*). (32) 

Now we will discuss the influence of massive fermions 
on a YM black hole. Yang-Mills equations of motion 
in the presence of fermions (after bosonization) are the 
following: 

V„[c2r']-2|w|2B'^ = AvV, (33) 

V„V«w~ 1) +2wB„B" = 

^■^sin(2^/^(^)). (34) 

From equation (l33T l we see that, even if we initially set 
Ba = 0, fermions give a nonzero contribution to the elec- 
tric part of YM field. This means that contrary to the 
massless case the presence of fermions leads to dyonic 
structure of a black hole. But as was shown in 01211 dy- 
onic Yang-Mills black hole necessarily has infinite mass. 
On the other hand, form equation (|34} we see that time 
dependent fermions give a nonzero contribution to the 
magnetic part of YM field. Through Einstein equations 
this leads to a time dependent line element and the de- 
struction of staticity of a black hole. 

CONCLUSIONS 

Now we shall present a short summary of our results 
concerning the influence of fermions on a primordial 
black hole modeled by a magnetic Yang-Mills black 
hole. 



First we will discuss the massless case. By asymptotic 
analysis we find an evidence that regular solutions to 
equations (fTST i and ( fT6] l should be time dependent. These 
equations describe bosonized massless fermions and give 
a nonzero contribution to a magnetic part of Yang-Mills 
field. Because their solutions are time dependent the re- 
sulting magnetic part of YM field will be time depen- 
dent and, through Einstein equations, elements of a met- 
ric tensor will also be time dependent. 

In massive case bosonized fermions are described by 
solutions to equation ( |29] l. The same type of analysis like 
in massless case also reveals the destruction of staticity 
of our black hole. Moreover, massive fermions will lead 
to the formation of a dyonic black hole. But, as was 
shown in [12], the presence of an electric part of Yang- 
Mills field leads to infinite mass of the resulting black 
hole. 

In conclusion we may say that the presence of 
fermions and their interaction with the considered PBH 
will lead to its destruction through mechanisms de- 
scribed above. 
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